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ABSTRACT 

In accretion systems, outflows may have significant influence on the luminosity fluctu- 
ations. In this paper, following the Lyubarskii's general scheme, we revisit the power 
spectral density of luminosity fluctuations by taking into account the role of outflows. 
Our analysis is based on the assumption that the coupling between the local outflow 
and inflow is weak on the accretion rate fluctuations. We find that, for the inflow mass 
accretion rate M (x r s , the power spectrum of flicker noise component will present a 
power-law distribution p(f) oc / _ ( 1 + 4;5 / 3 ) for advection-dominated flows. We also ob- 
tain descriptions of p(f) for both standard thin discs and neutrino-cooled discs, which 
show that the power-law index of a neutrino-cooled disc is generally larger than that 
of a photon-cooled disc. Furthermore, the obtained relationship between p(f) and s 
indicates the possibility of evaluating the strength of outflows by the power spectrum 
in X-ray binaries and gamma-ray bursts. In addition, we discuss the possible influence 
of the outflow-inflow coupling on our results. 
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1 INTRODUCTION 

The emission of Galactic Black Hole Binaries (BHBs) and 
active galactic nuclei (AGN) displays a significant aperiodic 
variability on a broad range of time-scales. The Power Spec- 
tral Density (PSD) of such variability is generally modeled 
with a power law, p(f) oc , where p(f) is the power 
at frequency /, and the power-law index j3 keeps a con- 
stant in a certain range of /, but changes among different 
ranges. At high frequencies, the PSDs of both BHBs and 
AGN present a steep slope with fi ~ 2. On the contrary, 
below a break frequency, typically at a few Hz for BHBs, 
they flatt en to a slope wit h /3 ~ 1, representing the flicker 
noise (see lKing et a!ll"2004l and references therein). 

Several models have been proposed in order to under- 
stand this nearly featureless c haracter of p ower spectra. The 
so-called "shot noise models" (|Terrelllll972l ) attempted to de- 
scribe the light curves as a series of independent overlapping 
shots with specific time-scales, amplitudes, and occurence 
rates. Due to lack of physical picture in this scenario, vari- 
ous physically motivated ideas have been put forward sub- 
sequently, such as the fluctuation s of hydrodynamic o r mag- 
netohydrodynam i c tur bulence (|Nowak fc Wagoned 1 19951 ; 
lHawlev fc Krolikl l200ll). magnetic flares or density fluctu- 
ation in the corona llGaleev et aflll995l; iPoutanen fc Fabian! 
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Lyub arskii's general scheme (|Lvubarskiil 1 19971 ; iKing et all 
12004) . In the Lyubarskii's scheme, it was noted that any 
variation of accretion rate, which is caused by small am- 
plitude variations in the viscosity, would induce a varia- 
tion in the accretion rate at the inner radius of the disc, 
where most of the energy is released. Moreover, observations 
showed that the variability is non-linear and the rms vari- 
ability is proportion al to the average flux over a wide range 
of time-scales (e.g.. Illttlev fc McHardvl l200ll ; lUttlev et all 
120051 : iGleissner et aJj|2004h . It indicates that the short time- 
scale variations are modulated by the longer time-scales, 
which favors the Lyubarskii's scheme. 

However, the observed power spectra often deviate from 
the form / . For example, the PSD of Cyg X-l is well de- 
scribed with the form / _1 in the soft state, however, ex- 
hibits the form /~ 1,3 in the hard state (|Gilfanovll2010t ). In 
particular, it is shown that the power-law index is around 
0.8 — 1.3 both in t he soft state of BHBs an d in narrow-line 
Seyfert 1 galaxies l|janiuk fc Czern"y||2007f ). Such a disper- 
sion of the power-law index reveals that there must exist 
some other mechanism. A radius-dependent amplitudes of a 
fluctuations may help to alleviate the discrepancy between 
theories and observations. However, it remains unclear why 
the fluctuations of a should have a strong radius-dependent 
form. 

In the present paper, we will take into account an- 
other mechanism, outflows, which is a popular phenomenon 
in accretion systems and has strong observational evidence. 
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One of the best examples comes from Sgr A*, whose center 
harbors a supermassive black hole surrounded by an accre- 
tion flow that is likely to be in t he form of the adve ction- 
dominated accretion flow (ADAF TNaravan fc Yilll994l ). Ra- 
dio polarization observations constrain the accretion rate in 
the innermost region is nearly two orders of magnitude lower 
than that measured at the Bondi radius fe.g.. lMarrone et al.l 
2006), which indicates that intense outflows may present in 
this system. Besides, the absorption lines from highly ionized 
elements, which have been detected in the X-ray spectrum of 
some microquasars suc h as GRO J 1655-40 (|Ueda et al.ll 19981 ; 
lYamaoka et £d] l200ll; fMiller et all l2006h . GRS 1915+105 
l|Kotani et al.ll2000l ; iLee et al.ll2002T) and A toll sources (see 
e.g. the review by iDi'az Trigo et all I2OO6I and references 
therei n) , also indicate th e existence of outflows. On the other 
hand, iJiao fc Wul |201ll ) found that outflows generally ex- 
ist in accretion discs no matter t hat the flow is advection - 
dominated such as the slim disc (|Abramowicz et aflll988l ) 
and the ADAF , or is radiation-dominated such as the stan- 
dard thin disc (jShakura fc Sunvaevlll973l ). In particular, for 
the three types of advection-dominated flows: ADAFs (gas 
internal energy dominant), slim discs (trapped photon en- 
ergy dominant), and hyper-accretion discs (trapped neu- 
trino energy dominant), outflows may be sign ificantly strong 
due t o positive Bern oulli parameters (e.g., INaravan et al.l 
1997; ILiu et al.ll201ll ) or the large radiation pressure (e.g., 
fc Lull2007l ). Furthermore, outflow s have generally been 
found in many simulation works (e.g.. lOhsuga fc Mineshigg 
l201ll and references therein). 

In the Lyubarskii's scheme, the power spectrum of lu- 
minosity fluctuations is sensitive to the varying mass accre- 
tion rate, thus we expect that outflows may have significant 
effects on the power-law index. The paper is organized as 
follows. In Section 2, we investigate the fluctuation power 
spectrum under a radius-dependent accretion rate following 
the method of Lyubarskii. The potential application of our 
results to observations is discussed in Section 3. 



2 LUMINOSITY FLUCTUATIONS WITH A 
RADIUS-DEPENDENT ACCRETION RATE 

2.1 Evolution equations of the fluctuations 

In the present work, following the Lyubarskii's general 
scheme, we revisit the power spectral density of luminos- 
ity fluctuations by considering a radius-dependent accretion 
rate. The relevant processes are the conservation of mass 
and angular momentum. With a radius-dependent accretion 
rate, these two processes are described as follows: 



9E _ J_ ( dM_ 

dt ~ 27rr ( dr ~ 



d(Amr 2 ) 
dr 



d 
- — (27IT T^) + $Z«, 
or 



(1) 
(2) 



where E, M , Q and T rv are the surface density, accretion 
rate, angular velocity and rip component of the stress tensor 
at the radius r, respectively. $ is the change of accretion rate 
over r in the stationary accretion state, i.e. <3> = dM /dr, and 
Z$ is the angular momentum of $. 



With Qr 2 = l in , F(= -T rv r 2 ) — at the inner radius 
ri n , we can obtain 

F = 7^- [A/fir 2 - Mn/in- / $l<s>dr] 



2tt' 

Afttr 2 
2n 



[1 - (A7i n /i n + / $l$dr)/(Mnr 2 )]. 



With the following definitions: 
h = fir 2 , 



(3) 



Q = 1 - (Minkn + / $Z*dr)/(Afflr 2 ), (4) 
the above description of F can be simplified as 



F 



Mh 



For standard thin discs and ADAFs, Q oc fi K = \/GM/r 3 , 
and thus h oc r 1//2 . We assume h = bV GMr in this paper, 
where & is a constant. If the accretion rate has weak de- 
pendence on the radius and oc Q,r 2 oc y/r, the value of 
Q would remain nearly constant. In order to simplify the 
problem, we assume Q to be constant in our analysis. The 
validity of this assumption is discussed in Appendix A. 

With x = (l/2?r) / $(Z$ - fir 2 )dr, equation © be- 
comes 

d[27v(F + x)} 



M 



dh 



(5) 



Substituting the above equation into the equation (1), we 
have 

9E _ b 4 (GM) 2 d 2 {F + X ) 
dt _ 



b GM 



<t>. 



(6) 



2h s dh 2 2nh 2 

The relationship between E and F can be deduced with a- 
prescription and the local balance between the heating and 
the cooling in the disc. The formula is similar to that with a 
constant accretion rate ijFilipov L.G.II1984I : lLvubarskii|[l987l ; 
INaravan fc Yilll994l ). i.e. 

b 4 (GM) 2 F 1 ~ m 

2(1 - m)Dh 3 -" ' ( > 

where the exponents m and n are determined by the disc 
model, which are discussed in §2.2, and D is a function of 
a. It should be noted that the analysis presented here is not 
strict for ADAFs. A stringent analysis for the ADAFs may 
refer to Appendix B. 

Substituting equation into equation ([6]), we have 



8F _ 2DF m 
~dt ~ b 4 {GM) 2 h n - 3 



,b A {GM) 2 d 2 (F + x ) 



2b? 



b GM 
2ixh 2 



$1 



dh 2 

F dD 
(1 - m)D~dt' 



(8) 



By assuming a = ao[l + P(t, r)], where /3(t,r)(<C 1) is de- 
correlated at different radial scales and its correlation time- 
scale is of the order of the local viscous time-scale, the dis- 
turbed quantities D and F are 



D = D {1 +77/3), F=( 



M h 
2ir 



where 77 is defined as D oc a v and the subscript denotes un- 
perturbed quantities. Substituting the above equations into 
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equation @, and including the following stationary condi- 
tion 

AGMf 9 2 ( ^ Q + X) _^£* = 0) 



2h 3 dh 2K 2tt 
we have, 



2irh 2 



chp _ D Q m M d 2 j, MphQ dy/3 

dt h n ~ m [ 2tt ' dh 2 2tt(1 - m) dt ' ( ' 

Here, we assume that the coupling between the local out- 
flow and inflow is weak on the accretion rate fluctuations. A 
simple discussion on this issue is presented in §2.4. 



2.2 The analysis of luminosity fluctuations 

Assuming the radius-dependent accretion rate is 



M = M ,i„( — y 

Tin 



(10) 



where Mo, in is the accretion rate at the inner radius n n , we 
have 



dt 



-,n — m — 2sm 



2tt 



dh 2 



Mo, infeQr? d h 2 s 



(11) 



2n(l~m) v 6 2 GMr in ' 

and C = D Q m /(GMr in b 2 Y m . This is a linear diffusion 
equation with a radius-dependent accretion rate. For s = 0, 
i.e. a constant accretion rate in the system, we have Q — 1 
and equation JTT} is reduced to 



di> Do ,M ,in. m d 2 i> , M ,inhrj d(3 



-( : 



+ 



(12) 



dt h n - m v 2tt ' dh 2 ' 2tt(1 - m) dt' 
which is the exact form of equation (9) in lLvubarskiil Jl997t) 



I n general, the solution of equ ation ()11[) is ( Lvubarskiil 
1 19971 , iLvnden- Bell fc Pringlelll974f ) 



47r(l-m) J 



t-V 



. (x + Xi)k K 2 XXl d — i xf ls 

CXP[ ' 4(t-t0 ]Il[ W^) ] dF m ' Xl \b 2 GMr in y 



where 



I = 



1 



2 + n — m — 2sm ' 



(13) 



k 2ty 
From equation (2), the accretion rate M(t,x) is 
M(t,x) - M(0,x) = 



_ x ±is 
G(t, x\ t', Xl ) — [p(t', xi) ,p GMrta y Wdx x , 



10 JO 

where 



G(t, x; t' , xi) — 



r)QK A x l - l x[ +1 



i 2 i 2\ 2 

[x + Xi)n 

Mt-v) 



1X2 eX Pi 



81(1 - m)(t - t') 

2 2 

{ r K XX\ . r Hi XX\ , , 

— ; r — xl i\—. r r, 
l 2(t-t') L 2(t — t') 



and I v (z) is the Bessel function of the imaginary argument. 
Based on the above equation, the power spectrum of M(t, x) 



can be obtaine d following a comp lex calculation as presented 
in Section 4 of iLvubarskiil (|l997h . Here, we directly present 
the result of the power spectrum and focus on the effects of 
outflows on the luminosity fluctuations. If \J < (ft) 2 > oc r^, 
the power spectrum p(f) of M(t,x) is 

P(f) oc /-ti-MW+'H, (14) 

which indicates the power-law index /3 = 1 + 4/(^ + s). In 
this work we ignore a fluctuations, i.e. £ = 0, and then the 
expression of /3 will be reduced to 



13 = 1+ Als. 



(15) 



For an advecti on-dominated flow, we have m = 0, 
rc = 1, and D oc a (|Naravan fc Yi|["l994l ') for equation (7), 
thus equation (13) indicates 1 = 1/3 and therefore ft can be 
simplified as 



(16) 



We would like to stress that the above formula should be 
valid for all the three types of advection-dominated flows 
mentioned in the first section. 

For a standard thin disc, it is well-known that there ex- 
ist three regions according to different dominant mechanisms 
for opacity and pressure. We have m = 0.3, n = 0.8, and 
D oc a ' 8 for the outer region (p ~ p gas , k ~ Kg), m = 0.4, 
n = 1.2, and D oc a ' 8 for the middle region (p ~ Pgas, 
k ~ Res)) and m = 2, n = 7, and D oc a for the inner region 
(p ~ Prad, K ~ K es ), thus we easily obtain the following from 
equations (13) and (15): 

40 



13 = 1 + 



25 - 6s 



for the outer region, 

P = l + ^s 

H 7 -2s 

for the middle region, and 
/3 = 1 + 



7 -4s 



(17) 



(18) 



(19) 



for the inner region. We would point out that outflows may 
be negligible in standard thin discs expect for the inner re- 
gion, which is radiation pressure dominated and may suffer 
the thermal instability. Outflows in the inner region may be 
significantly stronger than that in the other two regions due 
to the thermal instability. However, it remains under con- 
troversy whethe r the inner region is indeed thermally unsta- 
ble or not (e.g.. iHirose et al"T l2009'). and many mechanisms 
have recently been proposed to suppress the instability ( e.g., 
IZheng et al.fl201ll ; iLin et al.llioilt ICiesielski et~aHl201lh . 

Apart from the above mentioned photon-radiation- 
dominated flows, neutrino-dominated accretion flows 
(NDAFs) have a lso been w i dely studied (e.g. . 



Popham et aD Il999l: iNarayan et aD l200ll: Iciu et all 120061 : 
Chen fc Bcloborodovl 120071 : ILiu et alJ l2007h . which may 



acco unt for the cent r al eng ine of gamma-ray bursts (GRBs, 
e.g., iNaravan et all Il992f ). The outer region of NDAFs 
may be advection-dominated since neutrino cooling cannot 
balance the viscous heating due to low temperature and 
density. This re gion can be reg arded as an extension of 
slim discs (e.g., iLiu et aD l2008h . On the other hand, for 
high mass accretion rates such as m > IMq s _1 , the inner 
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cut radii (|Lvubarskiilll997l) . the power spectrum for the flow 
with M«r s can be written as 



p(/)«M 2 / -1 ocr 2 7"\ 



(21) 



where M is the accretion rate at the radius r which con- 
tributes to the flu ctuation at t he fre quency /, and the term 
/ _1 is taken from iLvubarskiil i|l997h for the constant accre- 
tion rate case. 

Taking ADAFs as an example, we have / oc l/t v i s and 
tvis = (r/H) 2 /(aQ) oc r 3/2 , i.e., r oc /~ 2/3 . Equation $2]} is 
therefore simplified as 

which is exactly the same form as equation (16). For other 
disc models, such as standard thin discs and NDAFs, the 
results can be understood in the same way. 



Figure 1. Variation of /3 with s for three types of accretion 
models: radiation-pressure-supported standard thin discs (the 
short dashed line), ADAFs (the solid line), and NDAFs (the 
long dashed line). The region of photon-cooled discs and that 
of neutrino-cooled discs are marked by the cross symbol and the 
dot symbol, respectively. 

region may also become advection domi nant due to the 
large optical depth for neutrinos (e.g., iDi Matteo et al.l 
12002] ; iGu et all 120061 ). Outflows may occur in NDAFs in 
particular for these two advection-do minated regions owing 
to positive Bernoulli parameters fe.g.. lLiu et al.ll201ll ). 

For NDAFs, with the balance between the cooling rate 
by pair capture and the ene rgy dissipation rate per unit 
volume (jPopham et al.l I"l999l . Eqs. 5.4-5.5), we have m = 
0, n = —0.4, and D oc a , thus equations (13) and (15) 
directly give the form of /3: 

P = 1 + \s. (20) 

The relationship between and s for all the mentioned 
accretion models are presented in Figure [T] where the short 
dashed, solid, and long dashed lines correspond to the in- 
ner region of standard thin discs, ADAFs, and NDAFs, re- 
spectively. The region with cross symbols represents photon- 
cooled discs (cooling dominated either by photon radiation 
or by advection of photon or gas energy), whereas the region 
with dot symbols represents neutrino-cooled discs (cooling 
dominated either by neutrino radiation or by advection of 
neutrino energy). The figure clearly indicates that, for com- 
parable s, a neutrino-cooled disc generally has a larger fi 
than a photon-cooled disc. In other words, j3 in a GRB sys- 
tem is generally larger than that in a BHB system. 

2.3 A simple analysis 

The dependence of /3 on s can be well understood with equa- 
tions (JS| and (|12|) . Following these two equations, the am- 
plitude of SM(r), which is the accretion rate variation pro- 
duced at a radius r, is proportional to the accretion rate at 
that radius, i.e., 5M(r) oc M(r). Since variations of the ac- 
cretion rate in the inner region may be represented as a sum 
of independent accretion rate variations produced at differ- 



2.4 Possible influence of the outflow-inflow 
coupling 

As mentioned in §2.1, our results are based on the assump- 
tion that the coupling between the local outflow and inflow is 
weak on the accretion rate fluctuations, i.e., the response of 
outflows to the inflow fluctuations is weak. In principle, the 
outflow should have fluctuations related to the variations of 
the local inflow accretion rate. However, the mechanism for 
outflows may be complicated, and therefore the fluctuations 
of the outflow may not be determined simply by that of the 
local inflow. For example, radiation from the inner region 
of a disc can heat up the materials i n the outer region and 
thus outflows can b e produced fe.g.. iBegelman et al.lll983l ; 
iMetzger et al.|[2008t ). In such case the outflow may be more 
relevant to the inner disc rather than the outer region where 
it occurs. 

Here, we make a simple discussion on the possible in- 
fluence of the local outflow-inflow coupling as follows. If the 
coupling exists, the amplitude of accretion rate fluctuation 
may vary while it propagates into the inner region. With 
the assumption that 8M(r)\ ri is the accretion rate varia- 
tion at the radius r\ induced by SM(r), we introduce a fac- 
tor D(r,r\) to generally describe the change of the above 
amplitude owing to the coupling effects, i.e., SM(r)\ ri — 
D(r,r\)8M(r). It is easy to find the form of D(r,r\) in the 
following two situations: [1], if the coupling is negligible, the 
amplitude of SM(r) will keep unchanged during its propa- 
gation, i.e., SM(r)\ ri = 5M(r), thus D(r,n) = 1; [2], if the 
local outflow is in strong coupling with the local inflow, the 
relative amplitude 8M (r)\ ri / M (r\) will keep unchanged for 
varying n, i.e., D(r,r\) = M(ri)/M(r). Obviously, our re- 
sults are under the former situation. On the contrary, for the 
latter one, the power spectrum of accretion rate variations 
at the inner radius ri n can be expressed as 

p(f) oc D^nnfM 2 /- 1 = Mir 1 oc r\ 

which is the sa me form as the re sult with a constant accre- 
tion rate, (e.g., Lyubarskii 1997). The detailed prescription 
of D(r, ri) is, however, beyond the scope of the present work 
due to the complexity of outflows. Nevertheless, we can ex- 
pect that a real flow may exist between situations [1] and 
[2], and therefore the value of j3 may be located between 
unity and the results presented in §2.2. 
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3 SUMMARY AND DISCUSSION 

In this paper, we evaluate the effects of outflows on the lu- 
minosity fluctuations with the Lyubarskii's general scheme. 
With a radius-dependent accretion rate M a r s , the power 
spectrum of the luminosity fluctuations is p(f) oc / , 
where the value of /3 varies with the disc structure. By 
assuming that the coupling between the local outflow and 
inflow is weak on the accretion rate fluctuations, we ob- 
tain the following explicit expressions of /3 for different 
disc models: j3 = 1 + 4s/3 for advection-dominated discs, 
/3 = 1 + 40s/ (25 — 6s) for the outer region of standard 
thin discs, /3 — 1 + 10s/(7 — 2s) for the middle region, 
ft — 1 + 4s/(7 — 4s) for the inner region, and f3 — 1 + 5s/2 
for NDAFs. The above expressions imply that in a GRB is 
generally larger than that in a BHB for comparable s. The 
expressions of j3 indicate the possibility of evaluating the 
strength of outflows by the power spectrum in X-ray bina- 
ries and GRBs. In addition, if the coupling is not negligible, 
the value of f3 will probably be located between unity and 
the value presented in the above expressions. 

In both BHBs and AGN, ADAFs are usually adopted 
to describe the quiescent state, the low/hard state, and the 
corona which lies above a cold disc. ADAFs may produce 
significant outflows, and therefore the power spectrum can 
deviate from /~ based on the present analysis. The exact 
value of s is, however, difficult to estimate from the theoreti- 
cal point of view, except for th e general constraint < s < 1 
l)Naravan fe McClintockll2008l). On the other h a nd, some ob- 
serva tions indicate s ~ 0.3 l|Yuan et al.l 120031 ; IZhang et alj 
|2010| ). Taking this value, we obtain (3 = 1.4 for ADAFs, 
which is close to 1.3, the power-law inde x of PSDs presented 
in th e low mass X-ray binary systems l|Gilfanov fe Arefievl 
2005). The quantitative difference may be relevant to the 
coupling between the outflow and inflow as discussed in §2.4. 

If there is only outflow that operates in the accreting 
system, s should be positive. However, s can also be nega- 
tive due to the evaporation mechanism of a cold disc. For 
the accreting black hole in BHBs, the observed power-law 
components in the X-ray spectra are generally attributed to 
hot, tenuous plasmas, namely accretion disc coronae. Due to 
the high temperature in the corona, the interaction between 
the disc and corona wo uld lead to mass e vaporating from 
the d isc to the corona (|Mever et al 1 l200d ; ISpruit fc Deufell 
2002). In this case, the value of s for the corona should be 
negative if outflows are not strong, and therefore it is quite 
possible for /3 to be less than unity. Consequently, in this 
scenario s for the underneath cold disc should be positive. 
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APPENDIX A: THE VALIDITY OF THE 
ASSUMPTION OF CONSTANT Q 

In this section, we analyze the validity of the assumption of 
constant Q. 

With a radius-dependent accretion rate described in 
equation |T0]), we have 



$ = s Mo4n ( _r_ )( .-i) 

fin Tin 

then equation ([4]) can be expressed as 
Q ~ 1 MVLr 2 

We assume that the specific angular momentum correspond- 
ing to $ is proportional to that of the gas in the disc, i.e. 



l<s> oc f2r oc r 



1/2 



This is the case for thermal energy drivi ng outflows, mag- 
netic field centrifugal accelerating wind (|Maver fc Pringlg 
l2006h and disc evaporation model. 

For the general situation (s 7^ —1/2), we have 



Q = 1 



M in l in + 



s+0.5 



ML 



MQr 2 



s + 0.5 Q.r 2 
and therefore 



( ' in S + 0.5'* 4n)( r^'nr 2 ' 



Q — > const, for s > —1/2, 

where l<s> t i n is the angular momentum of "I> at the inner radius 
ri n . Thus, the analysis presented in this paper holds for s > 
— 1/2. For s < —1/2, our analysis may present a qualitative 
result. 



9 , 1 (dM , . 

o 2 o 2 d P n 

pril - pril K ~ — = 0, 



2Hp 
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dt r dr 



r = XS- (2Hr 2 r rv ) 



2nr 2 



dr 



1 f dp dp p dp p dp 
72 — 1 V at dr p dt p dr 



where 71 and 72 are the usual generalized ratios of the spe- 
cific heat , r rv — —ap, H/r = const., and $, Al = Z$ — rv v 
maintain the value of the stationary state. The self-similar 
solution of the above equations is 

-1/2 -3/2 -5/2+s -3/2+s 

v r fl ocr ', SZq oc r , po oc r , po oc r ' . 

We introduce small deviations of the disc parameters from 
the stationary parameters as follows, 

a = a (l+P), P' = r 3 p, 

v T = — aictorflK (l + r~ s u) , = 02^1k (l + r~ a ujj , 
p — az*Jr£l%_r a (l + r~ s <5) , p = 04r~ 3 ^ 2+a (l + r~ a o} . 

where ai, 02, a-i, and 04 are constant. With the above equa- 
tions, the evolution equations of the perturbed variables are 

1 da d (v + a) 

— —air- 



dr 



aoO.K dt 

a^r — 038 — 2a4a\ui — 04 (a% — 1) a = 0, 

dr 2 y ' 

1 duj dui az 86 1 ,, „ . 1 

— — — - air— H r— - -ai (1 - 2s) w - -aiv- 

allK dt dr 0,40,2 dr 2 2 

1 a 3 03 d/3' a 3 -/ 
-aitr + S = r— + P , 

2 20402 0402 Or 20402 

1 dS 1 da dS da , . 

~FT "57 _ 7i-p^--77 - a ^ r 7T +7iair- 5 72 - 1) x 

ailK dt aiiK at dr dr 

5 3 \ . 3 
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dui 
a 2 r— + 
dr 



l<H + [»-2 7lOl + g (72 - 1) 0,2 

(h 3 \ / 3\ 

I - - s - -71 + S71 \ an> - s-y 1 a 1 a + I s + - 1 X 

3 

(72 — 1) a 2 u} = — - (72 — 1) a 2 /3'. 

The accretion rate is 

M = —4-Krv r pH, 

and its fluctuating component is 

m (r, /) = 47raia4Qo f — ^ [w + c] . 



S+ 



For W(/3) oc r°, the power spectrum is (see iLvubarskiil 
1 1991 Section 5) 

p(/)°c/- (1+ H 



APPENDIX B: FLICKER NOISE IN ADAFS 

The dynamic equations of ADAFs read as fo llows 
l|Naravan fc Yilll994l ; iKato et al.ll200sl ; lLi fc Caoll2009l ). 



